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Abstract 

Physical properties of anisotropic superconductors like the critical temper- 
ature and others depend sensitively on the electrón mean free path. The 
sensitivity to impurity scattering and the resulting anomahes are considered 
a characteristic feature of strongly anisotropic pairing. These anomahes are 
usuahy analyzed in terms of s-wave impurity scattering which leads to uni- 
versal pair breaking eífects depending on only two scattering parameters, the 
mean free path and the impurity cross section. We investígate here corrections 
coming from anisotropies in the scattering cross section, and find not only 
quantitative but also qualitative deviations from universal s-wave isotropic 
pairbreaking. The properties we study are the transition temperature, the 
density of states, quasiparticle bound states at impurities, and pinning of flux 
lines by impurities. 

PACS. 74.20.-z, 74.20.Fg 
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I. INTRODUCTION 



Unconventional anisotropic pairing is evidently realized in high temperature supercon- 



pairing is an important tool in analyzing the symmetry of the pairing amplitude, and is the 
subject of a number of experimental and theoretical wor Most of calculations were done 
assuming an s-wave impurity scattering potential Uimp{p, p') = const, and taking either the 
Born limit {N{0)uimp <^ 1) or unitarity limit {N{0)uimp ^ 1). 

Surprisingly, a number of experiments on the optimally HTS oxides have shown that d — 
wave pairing is quite robust, i.e. not very sensitive to various kinds of impurities and defects. 
For instance, the decrease of the critical temperature Tc{pimp), with increasing residual 
resistivity, Pimp, is much smaller than the theory with the s — wave impurity scattering 



effects by impurities in HTS oxides was proposed by the authors of Refs.El,H, who invoked a 
momentum dependent impurity scattering potential with an appreciable contribution in the 
d-channel. The microscopic theory in Ref.i accounts for the renormalization of the impurity 
potential by strong correlations, which gives rise to a pronounced forward scattering peak, 
while backward scattering is suppressed, as first proposed in Ref.S (see alsJ). Application 
of this theory to impurity scatteringi shows that in addition to the contribution in the 
s-channel there is a significant contribution to the Born amplitude from the d-channel of 
the same magnitude , in particular for low (hole) doping concentration, 6 < 0.2. As a 
consequence, the decrease of Tc{pimp) with increasing p¿mp is much slower than the theory 
with exclusively s-wave impurity scattering predictíl. This renormalization effect explains 
the robustness of d-wave pairing in HTS oxides. 

One may rise the question whether this robustness also holds far away from T^ and for 
very strong scattering potential, for instance in the unitarity limit. To answer this question 
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we shall analyze a class of models by calculating the scattering T-matrix with an impurity 
potential that depends on the scattering angles. 

A related class of problems, which we study in section III, is related to the impurity 
scattering in the two-band model. Recently, several models for the pairing mechanism in 
HTS oxides based on two-band and multi-band modelsi,0,0,lli were suggested, and impurity 
effects studied in Born approximation. Magnetic and non-magnetic interband scattering can 
lead in this model to a lowering of the critical temperature and also to a relative sign change 
of the order parameters in different bandJi^. In section III we analyze the changes in the 
two-band model when going beyond the Born limit. It was shown that in the unitarity limit 
the Anderson theorem holds. 

In previous sections we studied a homogeneous superconductor with homogeneously dis- 
tributed impurities. Selected inhomogeneous problems are studied in section IV, such as the 
bound states at an impurity, and the pinning energy at an impurity (defect) of singly- and 
double-quantized vórtices. 

II. ANISOTROPIC SCATTERING IN ANISOTROPIC AND HOMOGENEOUS 

SUPERCONDUCTORS 

In the following we analyze superconducting properties of anisotropic superconductors in 
the presence of momentum-dependent nonmagnetic impurity scattering by the quasiclassical 
equations of Eilenberger, Larkin-OvchinnikovllÍ,lii {ELO equations ). For a homogeneous 
distribution the quasiclassical Green's function matrix, 5f(pi?, R, cu^), is independent of R, 
and the quasiclassical equations read 




(1) 




(2) 



We assume weak-coupling superconductivity with Á(pj7)(= '¿A(pp')f2), where A(pp') is 



real. The 2x2 matrices tq = 1 and ti^2,3 are Nambu-Gor'kov matrices. The effect of 
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nonmagnetic impurities is described by the self-energy, ájmp, given in terms of the forward 
scattering part, í(pj?, p^, a;„), of the T-matriyl^ 1^ 



'5'¿mp(PF,í^n) = CÍ(pF,PF,í^n), (3) 

where c(<C 1) is the impurity concentration. 

For simphcity we assume an isotropic Fermi surface but pairing and impurity scattering 
are angle-dependent, i.e. A{pf) = A(s), ^(pi?,u;„) = g{s,n) and t{pF,p'p,LJn) = t{s,s',n) 
where s = pp/pp.The T-matrix is the solution of the equation 

í(s, s',n) = u{s, s')Í + A^(0) J ds"u{s, s")g{s", n)í(s", s',n), (4) 

where for the 2D systems, which we consider here, one has / ds{..} = J {..}d6/2'K. Since 
A(s) is real one has g = g2f2 + gsfs and i is given by t = tofo + + hT'2 + tsT's- 
Because the unperturbed solution has the form {un = T^T{2n + 1)) 



í/(s,^) = 1 ^ (6) 



then g{s,n) is searched for in the form 
where 

a)„(s) = Unis) + icits^s, s,n) (7) 

Á(s,u;„) = A(s) - ic¿Í2(s, s,n). (8) 

The self-consistency equation for A(s) is given by 

A(s) = N{0)TY: Í ds'V{s,s')g2{s',n), (9) 

where the pairing potential V^(s,s') =VpY{s)Y{s') is assumed in the factorized form with 
< Y'^{s) >s= 1. The latter implies that the order parameter has the form A(s) = A ■ Y{s). 



For convenience we define r„(= 07^) = c(7rA^(0))~^ and v{s,s') = 7rN{0)u{s,s'). In what 
foUows we consider the effects of anisotropic impurity scattering on the anisotropic pairing 
where < Y{s) >s= 0. 

1. Anisotropic impurity scattering and nodeless anisotropic pairing 
First, we consider the nodeless d-wave like pairing A(s) = A-Y{s) which is characterized 
by (y(s))s = O and l^^(s) = 1. This means that there is a finite gap everywhere on the 
Fermi surface, i.e. A(s) 7^ 0. It is interesting to mention, that besides the simphcity of this 
kind of pairing and its adequacy in some quahtative understanding of d-wave pairing it also 
appears to be a solution of the spin-bag mo dellli for HTS oxides. In this model the nodeless 
d — wave like pairing is due to residual (longitudinal and transverse) spin fluctuations on the 
antiferromagnetic background, where the AF order is distorted locally by hole doping and 
the spin-bag is formed around doped holes. The impurity scattering potential is assumed to 
have the form 

v{s,s') = vo + V2Y{s)Yis'), (10) 

i.e. it contains an anisotropic contribution in the same channel as the unconventional pairing. 
The solution of -E5'.(§) for Í3 and ¿2 is searched in the form 

Í3(s, s') = [Í3o(n) + Í32(n)F(s)r(s')]^3, (11) 

t2{s,s') = Í2{n)[g2{s,n) + g2{s',n)]. (12) 

(Note, in this model g2{s,n) = g2{n)Y{s), g2{s,n) = gliji), g^{s,n) = g^iji) and due to 
i?g.(^ one has g1{s^n) + g2{s,n) = —1.) The solution is given by 

i 

'{l + vi){l + vi) + {vo-V2rg'2{ny 



1 + V 

Í3o{n) = luvl—^ — — 2rTT \2~2í \ ^ (13) 



= ^'"•'°'Nl + «J)(l + ,.i) + (.„-..)^9Í(n)' 
while Í32(^) = t3o(^, t'o ^ ^2)- Several interesting results comes out in this case. 



(a) The critical temperature Te 

In the limit T —>■ Te Eqs. (y, |1^ — 11^ become 



where the pair-breaking parameter Fp^ is given by (r„ = C7r/A^(0)) 



(16) 



Note that T^ vanishes for F^^ 0.88Tco and this pairing is in some respects similar to d- 
wave pairing. It is apparent from -Eq'g. (p!5| — p!6|) that the pair-breaking effect of impurities is 
weakened in the presence of momentum-dependent scattering and it is even zero for = V2- 
The latter resuh has been previously derived in the Born approximationü. For f 2 ~ vq the 
slope dTc/dpimp can be very small even for appreciable valúes of pimp ~ F^^ = F„((To + d'2), 
because in that case Tpb -C Fí,, as indícate the experimental results oíi. The parameters ai 
are given by 



V? 



l + V. 



2 ' 



0,1,2... (17) 



The resistivity, p^mp, and the reduction of Te due to impurity scattering, Tc{pimp), depend 
on the classical transition rate, W{s,s') = F„ | íjv(s, s', n) pin the normal state0. This 
transition rate comprises all the needed Information on impurity scattering for either solving 
the normal state Boltzmann equation to determine pimp or the linearized gap equation Eq.(^ 
to determine T^. For the latter purpose one needs linear (integral) equation for (72 (s, n) which 
reads 

I cj„ I ^2(s,n) - A(s) + J ds'W{s,s')[g2{s,n) - g2{s',n)] = O, (18) 
where the normal state t-matrix ti\i{s,s',n) is the solution of the equation 

t^is, s' ,n) = v{s,s') — ÍTrsign{uJn) J ds"v{s,s")tN{s" ,s' ,n). (19) 
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Henee, measurements of Tc{pimp) curve carry not enough informations on the microscopio 
scattering data, i.e. the scattering T-matrix. More such informations are contained in spec- 
troscopic data on anisotropic superconductors at temperatures T <^ Te, such as tunnehng 
data or optical data at about gap frequency. 

(b) The density of states, Niu) = A^(0) Im / ds g^{s,iujn ^ u — ir]), depends in the 
presence of pair-breaking impurities significantly on the valúes vq and vi. It is knownlil that 
in the case of s-wave scattering only {v2 = 0) one has N{uj = 0) 7^ O for TuVq > A, and 
the highest valué, N{tü = 0) = A^(0)/[0.5 + 0.5(1 + (2A/r)2) V2]i/2, where T = T^ao, is 
reached in the unitarity limit. On the other hand, one obtains in the limiting case, vq = V2, 
a restoration of the gap, N{uj = 0) = 0. Despite of the strong scattering limit N{uj) is 
BCS-like. 

2. Isotropic impurity scattering and d — wave pairing 

Let US study a two- dimensional superconductor with the pairing function A(s) = A{(p){= 
A • 1^2 (v)) ~ cos2(y9 - d — wave pairing. Note this case seems to be more realistic for HTS 
oxides than the previous one, because the "eos 2(f" pairing has nodes at the simply connected 
Fermi surface. We assume that the isotropic impurity potential depends on the transferred 
scattering angle 

v{ip, ip') = Vq + 2vi COs((y9 — (f') + 2t>2 COS 2{ip — (f'), (20) 

where v{(p,(p') contains the pairing channel (~ 1^2 (v^) ^2 (v^')) ^oo- (This problem but with 
f2 = O is studied inEl but there is an inappropriate sign in the Í2-Hiatrix, which in fact 



corresponds to a magnetic impurity scattering). From £'gs.(|T^ — [T^ one obtains the pair- 
breaking parameter Vpb 

(l_^)2 + ^2(i^ 2) 

= r.[ao + ai], (21) 

where a = V2/V0 and a^, are given by -Eg.([r^). 

In order to analyze Tc{pimp) dependence, where the residual impurity resistivity pimp ~ 
Ttr, we need the transport scattering Ttr which is in this case given by 



Ttr = Tuial + 2cT¡ + 2a¡ - 2íTo[ai(l + ^) + ^2(1 + ^)]}. (22) 

VqVi VqV2 

If one wants to interpret depairing effects of impurities and robustness of pairing in HTS 
oxides in terms of the above results then the experimentsi imply that the ratio, Tpb/Ttr, 
should be minimum. In the case when V2 <C vi one obtains, Tpb/Ttr = 2, in both, the 
Born and unitarity, Vq,Vi — » 00, hmits. For, Vi <^ V2, the pair-breaking parameter, Tpf,, is 
minimized for a = 1/2 which gives, Tpf,/Ttr ~ 1/3, in both hmits. This means that the latter 
case is more appropriate candidate, than the case, V2 <^ fi, for the quahtative explanation 
of robustness of d-wave pairing in HTS oxides. 

III. TWO-BAND MODEL WITH NONMAGNETIC IMPURITIES 

The interest in two(multi)-band models and in the impurity effects is renewed after the 
discovery of HTS oxidell,0,0,0, where various kinds of the intra- and inter-band pairing 
and impurity scattering are considered. By assuming only intra-band pairing (a = 1, 2) 
the effect of the nonmagnetic impurities in the Born approximation is described by the 
equations (n enumerates Matsubara frequencies) 

Á„„ = A« + ^ ^lap (24) 

p /3n 

A« = 7rT Y. V=^, (25) 



I3,n 



in 



where Qan = \J^an + ^an' 7a/3 = ^a(3-^Á^) noumaguetic impurity scattering and 

^ap = ^1/3^/3(0) corresponding couphng constants. In Ref.0 are considered various 
possibilities for the suppression of the critical temperature, as well as the relative sign of 
Al and A2, in the Born limit for nonzero valúes of \af3 and 70,/?. We note some interest- 
ing conclusions obtained in Born approximation which shall be compared with the results 
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obtained in the unitarity limit: (i) the diagonal scattering rate 711 and 722 disappear from 
the linearized for T^, {ii) in the case An 7^ O, A22 = A12 = A21 = O the depression 

of Te is given by 5Tc/Tc = — 7r7i2/8Tc; (ii) for An = A22 7^ O and A12 = A21 = A_l < O one 
has sign{Ai/ A2) = —1 and 6Tc/Tc = —71(712 + 72i)/8Tc, while the sign of Ai and A2 is 
unchanged by the impurities. 

The í-matrix equation in the two-band model has the form {a, /5, 7 = 1, 2) ( we consider 
a rather small impurity concentration and neglect an interband hybridization) 

t(n) = ú + ^úN(0)g(n)t(n), (26) 

7 

where t{n) = I]^=o*í®^í' ~ g3®^3 + g2®^2 and (g) is the direct product of matrices in 
the band space (bold) and in the Nambu space (hat). g2, gs and N(0) are diagonal matrices 
in the band space. 

In the case of nonmagnetic impurities one has = u^®to and since gi = O one has 
<(n) = + u^N(0)[g3(n)tf (n) + g;{n)t^ (n)] 



tf (n) = u^N(0)Hg3(n)t^(n) +zg2(n)tf (n)] 
t^(n) = u^N(0)fe3(n)tf (n) +g2(n)t^(n)] 

t^{n) = u^N(0)[g3(nX(n) -zg2(n)tf (n)] (27) 

Let US consider for simplicity the case when = O but interband scattering, = 

= u O, and introduce three parameters 

^ 7r^iVi(0)iV2(0)^^ 



and 
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After some straightforward calculations one obtains the renormalized frequencies, üin, and 
order parameters Aj„ 



UJin = UJn + Ti T— ' (30) 

det 1 



ZAin, — ZAi + i 1 — — , [óí) 

det 1 



where 



det 1 = 2((T - l)crV'^ln + Af„(Ai„A2n + UJinU2n) 



[2{a - l)a + l]{Col + Kn)^ ^In + AL- (32) 



The solution for the second band is obtained from i?gs.(|3^ — ^2]) by replacing 1 <í=^ 2. In 
the Born hmit one gets 

UJln = UJn + i-l I , (33) 



+ a: 



2n 



Ái„ = Al + Ti ^^^^ , (34) 

V'^2n + AL 

i.e. the interband scattering mixes both bands. In the unitarity hmit a ^ \ {u ^ oo) the 
bands are decoupled, i.e. 

^an = í^n + Ta — ; „ = (35) 

an ' an 

Á„„ = A„„ + Aan ^^ ^gg^ 

V'^L + AL 

So, in this case the Anderson theorem is restored, i,e, the thermodynamic properties are 
impurity independent. 

The latter resuh can be generahze to the case 
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u 



N 



( \ 

au u 



V 



(37) 



u U22 

For u ^ 00 but a and U22 finite, a and U22 drop out from equations and the bands are 
decoupled with ujan and Aq,„ given by -Egs.(p5| — 133). At Te one has 



^In = + TiSigniUn) 



(38) 



^In I I ^2n 

From i?g.(^) it is seen that in the unitarity hmit, a — > 1, the renormahzed order parameters 
are decoupled and is unrenormahzed. For a < 1 it can be easily shown that is reduced. 



IV. SMALL ANISOTROPIC DEFECT IN ANISOTROPIC SUPERCONDUCTORS 

In what follows we consider the effect of a single impurity (small defect) with small scat- 
tering length a, which is supposed to be much smaller than the superconducting coherence 
length, I a |<^ .^o- Henee, the impurity can be considered as a localized perturbation, but 
with negligible renormalization of Á(pj7,R), giving rise to the quasiclassic equationS,Íll 

[{iUn + e\F ■ A(R))f3 - k{^F, R.), (5^(PF, R, í^n)] + Í^F^H^giVF, R, C^n) = 

= [í(pF,PF,í^n),5'¿mí(PF,R.,í^n)]<5(R.-R.¿mp)- (39) 

Here, 5^(pir, R, ciJ„) = ^(pi?, R, cUn) — ^jmí(PF, R, í^n)- The extra term proportional to 
5(R — Rfmp) describes a jump in ^(pi?, R, c<j„) at the site Rfmp of the impurity (defect), 
while the intermedíate Green's function ^¿^¿(pp', R, cu^) describes the quasiclassic motion in 
the absence of impurity (defect) and it is the solution of -Eg.(p^) by putting the right-side 
to zero. ^¿mí(PF, R, i^n) is normalized according to -Eg.(0). The t-matrix is the solution of 
i?g.(^ where 5'(pf, ^n) is replaced by gimt{PF, R = Rimp, ^n)- The change of the supercon- 
ducting free-energy in the presence of a single impurity (defect) is given by@,Í3 
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(5F(R,„p) = NiO)TY^ l^dxj^l d'RTr[5gÍPF,R,u^)A,ÍpF,R)], (40) 

where Ab(pir,R) and the vector potential A¿(R) are calculated in the absence of the im- 
purity. The Green's function, óg^ppjUjUn), must be evaluated for an order parameter 
Á(p^,R) = AÁb(p^,R). 

In the following we study the consequences of anisotropic impurity scattering for three 
selected examples of inhomogeneous anisotropic superconductors. 

1. Bound states due to the anisotropic impurity 

Let US consider the local change of superconductivity in the presence of a single 
anisotropic impurity with the potential v{s,s') given by -E'g.(0) and analyze the impurity- 
induced quasiparticle bound state and the change in the free-energy (5F(Rjmp)- By assuming 
that 27rd"¿ ^ Ep / Aq , where 2 = 0,2 and ai = vf / {1 + vf), the t-matrix is given by the same 
expression as -Egs.(|l3| — 0), but with g2{'n) is replaced by g^\n). The bound state energy 
^B,anis < Aq, which ís due to the pair-breaking impurity effects, can be obtained as a pole 
of the í-matrix which gives 



^B,anis — ^Oy^ — O'pb, (41) 

where 

- - - - ^2)^ . , X 

^0^2 

In the unitarity limit for both channels, i.e. fo ^ 1, f2 ^ 1 but V2/VQ finite, one has 
^B,anis coutrary to the unitarity limit for the s-wave scattering (fo ^ 1, 1^2 = 0) where 

^B,iso 0. However, the zero-energy bound state ujB,anis O appears when V0V2 = —1, i.e. 
if one channel is in the unitarity limit the other one must be in the Born limit. 

Due to the bound state there is a change (increase) of the free-energy 5F(Rjmp) = íi^imp- 
By solving Eq.{^% with gim,t{pF,R,^n) given by -^(/.(^ and t given by i^g^.j!^ — |14D one 
gets óFimp from Eq.{^ 

OTimp - J- 2^ UAOpb^ ^ \2A21L,2 _L,.,2 



12 



= 2Tln 



cosh(Ao/2T) 



(43) 



COsh[(l - Orp,)l/2Ao/2r] 



where apb is given in -Eg.p^). It is seen that there is a loss in the condensation 
energy,5F(Rjmp) > O, which is related to the pair-breaking effect of impurity. For vq = V2 
such an impurity does not affect superconductivity and 5F(Rjmp) = 0. 

The obtained results tell us that in for angle-dependent impurity scattering even a strong 
impurity potential may have very weak effect on Te, the bound state, and the free-energy of 
anisotropic and unconventional pairing. In that case the anisotropic pairing is robust in the 
presence of impurities. 

2. Pinning of single- vortex by a small anisotropic defect 

Because in HTS oxides strong correlations give rise to strong momentum-dependent 
charge scattering processes it is interesting to analyze the elementar y-flux-pinning potential 
of a small defect by using the approach of Thuneberg et al., , who showed that in s-wave 
superconductors the pinning energy of a small defect (a <^ ^o) is dominated by scattering 
processes at the defect. It is proportional to the product of the scattering cross section and 
coherence length (oc a^^o); instead of (naively believed) a^. The case of anisotropic pairing 
with s-wave impurities and near Te was recently studied inü. 

In what follows we study the effect of scattering anisotropy on the pinning energy of 
a small defect in an anisotropic superconductor at any temperature below T^. We use the 
model potential given in -Eg'.([To|) and assume that the vortex is placed at the defect. In order 
to calcúlate the elementary flux-pinning energy one has to solve the quasiclassical equations 
for various ballistic trajectories with R-dependent the vector potential, A(R), and order 
parameter. 



In the gauge where 6 is the angle with respect to the X-axis then A(R) has no radial 



For a qualitative discussion we will adopt a simplified vortex moderna which neglects the 



(44) 
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suppression of the order parameter in the vortex core and sets | Afc(pi;',R) |= Ao(pf), i-e. 
independent of R. Henee, the order parameter along a trajectory passing through the vortex 
center has constant magnitude but its phase changes abruptly by tt when going through 
the vortex core. This "zero-core model" gives the right order of magnitude of the pinning 
energy, 6Fpin(Rimp), when compared with the numerical calculations0. In order to calcúlate 
6Fpin{Rimp) two quantities are needed: gimt{PF,'R = Rimp,uJn) = gvÍPF,'R = Rimp,uJn) in 
the presence of the zero-core vortex and the impurity í-matrix calculated with the Green's 
function, ^„(pi7,R = Ilimp,uJn), of the zero-core model. The solution is straightforwardill 
and gives 

^„(pi7',R = Rimp,uJn) = —[{-A2T1 + AiT2)Y{9) + (-íia„)'r3], (45) 

and 

Í{Pf, Pf, UJn) = = -i'-iuOín^ni o »2 + ^~^~^]^3- (46) 



Here, = yct;^ + Ag. -5'q'.([39|) can be solved by the Fourier (or Laplace) transform which 
gives the expression for the pinning free-energy 

í^P» = íi^i-f ^^(^0, a,) + 5F^^{ap,) (47) 



5 = -2T In ^osh(Ao/2T) 

cosh[(l-a,,)V2Ao/2r]- ^^^^ 

-Egs.(^ — imply that, 5Fpj„ < O, and the vortex is attracted (pinned) by the defect. A 
comparison of Eq.(^l\) with the corresponding results for s — wave superconductors with an 
s-wave scattering potential shows, that in the former case two additional terms are present. 
The first one, depending on d"2, appears also in s-wave superconductors with anisotropic 
scattering accounted for. In fact óF^l^^^^ describes the reduction of the superconducting 
stiffness in the presence of impurities. For instance near i?g.(^) gives 
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SF^"' = -(-0 + • ^ - -7-6^^0^ ■ TAUn (50) 



and 6Fjj¡l^^^^ is proportional to the total scattering amplitude co + (72. For vortex far away 



pm 



from the impurity there is loss in the condensation energy 6F^^ [aph) due to pair-breaking 



effect of the impurity, i.e. 5F^^{apb) = —SFimp{crpb) where SFimpicpb) is given by Eq. 



Therefore this part enters in -Eq'. (|47|) with the negative sign, thus increasing the pinning 
energy when vortex is sitting on the defect and stabihzing it additionally. Near one has 



SFl^t\^eff) = -^p,-^, (51) 



For the s-wave scattering only, V2 = O, one has a2 = O, apb = ctq, and the pair-breaking effect 
is maximal while the condensation energy is gained maximally for vortex sitting on the 
defect. However, for, V2 = vq the pair-breaking of impurity is absent apb = O and 6F¡,t^ = O, 
i.e. in this case the pinning by the small defect is similar to that in s-wave superconductors. 

The physical picture of the vortex pinning by small defect given above is based on 
the known results based on the microscopic derivation of the Ginzburg-Landau equations 
in the presence of impurities. An explanation based on the quasiclassical approach is 
given in0,Íil and we briefly discuss it in order to develop an intuition for the case of a 
double-vortex pinning, which is studied below. Because the order parameter changes its 
phase by vr along the trajectories across the vortex core it leads to the phase change of 
gv{PF, R-, ^n){= dimtiPF, R-, ^n)) oji thc distaucc ^q, thus causiug a cost in the condensation 
energy, i.e. the maximal increase of the free-energy. Note, the function ^„(pj?, R, a;„) de- 
scribes the quasiclassical motion of particles (or pairs) along trajectories across the vortex 
core where the maximal phase change (vr) occurs. In the presence of defect the motion 
of particles is described by the function 5í(pj?, R, a;„) which contains scattering of particles 
to new directions where the phase change (mismatch) is less than vr and it costs less con- 
densation energy. Therefore the vortex is attracted to the defect because scattering helps 
superconductivity to sustain abrupt changes in the order parameter. The latter explains the 
contribution 6f!^¡I^'^^^ , while in the anisotropic superconductors, due to pair-breaking effects 
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of impurities, there is a gain of condensation energy —SFimp{creff) for vortex sitting on the 
defect. 

3. Pinning of double-vortex by small defect 

We extend the calculations to the pinning of multiply-quantized vórtices on small 
defects. First, a s-wave superconductor is considered and we put the question - is it possible 
to pin the double-flux-vortex ($ = 2$o) by the small defect, which is for simplicity charac- 
terized by the parameter ctq for s-wave scattering only? The "zero-core model" is assumed 
again. In that case the order parameter can be parametrized in the form 

Afe(p^,R)=| A(p^,R) le^*^ (52) 

For particle motion across the double-vortex core the order parameter does not change phase 
and in that case the solutions for gimtiPr, R- = R-imp, ^n){= g2v{PF, R- = R¿mp, í^n)) and for 
ÍÍPf, PF,uJn) are given by 

'i 

g2v{PF, R = R-ímp, ^n) = — [Aiñ + Asfs + (-í^„)r3] (53) 

and 

ÍÍPF, PF, í^n) = h% = -iUnanlu o i^- a 2 ^3, (54) 

where = (Tq/vq. Then by solving Eq.(^^ and by using -Eg.(^OD and this solutions one 
obtains the pinning energy of the double-vortex within the zero-core model 

«.,„^2rE/>A-^^^gl^>0. (55) 

The main conclusión coming out from -Eg.([55|) is that because SFpin > O the double-vortex 
in s-wave superconductors is repelled from the defect - i.e. the zero-core double-vortex can 
not be pinned. Contrary to the single-vortex, where the defect scatters particles to new 
directions where the phase change is smaller, in the case of double-vortex the particles are 
scattered to directions where the phase change is larger. However, it might be that the above 
obtained results are an artefact of the "zero-core model", where there is no suppression of 
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the superconducting order due to the vortex core, and numerical calculations are required 
for a realistic double-vortex structureü. 

In the case of an unconventional pairing, like that in Section J/.l, the order parameter 
is given by Afc(pi?,R) =| A(pj7,R) | exp{2i9)Y{9) and the t-matrix contains also terms 
^1,^2 7^ O leading to a decrease of the jump [i{pF,pF,^n), gimtÍPF,'^,^n)] in -Eg-P^)- In 
this case the pinning energy contains the additional term (gain in energy) due to the pair- 
breaking effect of the impurity, —SFimp, i.e. óFpm = 5F2v,pin — ^Fimp- Since 5F2v,pin is less 
positive (repulsive) than for a s-wave superconductor in -Eq'.(|55|), and because —SFimp < O 
one can happen that óFpm < O and even the zero-core double-vortex can be pinned by 
the defect. A realistic calculation of óFpin for anisotropic superconductors with anisotropic 
scattering of single and double-vortex will be discussed elsewhereü. 

In conclusión the anisotropic impurity scattering gives rise to new qualitative effects in 
unconventional and anisotropic superconductors, where for instance it " screens" the strength 
of the scattering in some quantities (like T^. - robustness of pairing, bound states, pinning, 
etc.) even in the unitarity limit. It seems that this situation is partly realized in HTS oxides 
where the d-wave pairing is robust in the presence of even very strong impurity scattering. 
In two-band models nonmagnetic impurities do not affect thermodynamic properties of s- 
wave superconductors in the unitarity limit for the interband scattering, contrary to the 
Born limit, i.e. in this case the Anderson theorem is restored in the unitarity limit. 
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